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Tjj- \ inating in the — K mass matrix, was discovered nearly forty years ago 



Kaon decays continue to provide invaluable information about the ap- 
proximate discrete symmetries of nature. CP- violation in ^ vrvrlQ, orig- 



^ ' and direct CP- violation in K-decays has been unambiguously established 0, 

^ ■ 0,0], through a recent remeasurement of Re (e'/e) = (28.0±3.0±2.8) x 10~^ 

O . by the KTeV collaboration 0. KTeV has also observed and studied^ the 

! rare decay — >■ Ti^TT~e^e~ . A large CP- violating asymmetry, i?cp = 

13.6±2.5±1.2 %, constructed from the final state particles was measuredQ, 
P^l consistent with theoretical predictions]^, ^ ^ This decay is dominated 

p by a one-photon intermediate state, — > 7r+7r^7* Tr+vr^e+e^ and Bqv 

receives a sizable strong interaction enhancement. 

A long standing problem in better understanding K-decays and a road- 
block to more precisely constraining the standard model of electroweak 
^ ' interactions or uncovering new physics is our present inability to compute 

d • the hadronic matrix elements of most electroweak operators to high preci- 

sion. The lattice provides the only direct method with which to determine 
these matrix elements, however, it is presently far from being able to com- 
pute matrix elements between multi-hadronic initial and final states. Chiral 
perturbation theory, xPT, is a framework in which the low-energy strong 
interactions of the lowest-lying pseudo-Goldstone bosons can be treated in 
perturbation theory. The external momentum and the light quark mass 
matrix are treated as small expansion parameters when normalized to the 
chiral symmetry breaking scale, A^ ~ 1 GeV. This article presents the 
xPT analysis of K^^ 7r+7r~e"'"e~, focusing entirely on the one-photon 
intermediate state, as shown in fig. (|I]). 

The matrix element for 7r'*"7r~e+e~, assuming CPT-invariance, is 

*Talk presented at the Kaon99 Meeting, University of Chicago, June 1999. 
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Figure 1: The one-photon intermediate state dominates 
The solid circle denotes the —>■ t:~^t:~j* vertex. 
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written in terms of three form factors G, F+ and F_, 
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where A;+__ are the positron and electron momenta respectively, q — k^-\-k_ 
is the photon momentum and p+ _ are the 7r+'~ momenta respectively. 

is Fermi's coupling constant, Si is the sine of the Cabbibo angle, / 
is the pion decay constant, and a is the electromagnetic fine structure 
constant. The form factors are functions of hadronic kinematic invariants, 
e.g. = F^(q^, q ■ p+, q ■ The smallness of Re (e'/e) suggests that to 
a very good approximation direct CP violation that may contribute to this 
decay can be neglected. For our purposes the only CP violation that will 
enter into this decay is due to e, indirect CP violation introduced by the 
wavefunction. In terms of the eigenstates of CP, Ki^2: the wavefunction 
is 



\Kl) = \K2) + e\K^) 



K 



, (2) 



where GP\K^) = +\Ki), GP\K2) - -IK^) and e = 0.0023 e*^^ . As 
direct CP violation is being neglected it is convenient to determine the 
contributions to G, F+ and F_ from Ki and K2 independently as the two 
contributions do not interfere in the total decay rate, F, or differential decay 
rate dV/dq^. The CP-odd component of the wavefunction, K2, gives 
contributions to the form factors with symmetry properties G — > +G, and 
F± under interchange p± ^ Pzp, while the contributions from the 

CP-even component of the wavefunction, Ki, have symmetry properties 
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Figure 2: The leading order contribution to 7r+7r~7* in xPT . A 

solid square denotes a weak interaction. Only the Ki component of the 
wavefunction can contribute at tree-level and this contribution is suppressed 
by a factor of e. 



F± under interchange p± — >■ (where it is understood that the 

interchange p± p^p also occurs for the arguments of the form factors). 

The lagrange density that describes the leading order strong and As = 1 
weak interactions of the lowest-lying octet of pseudo-Goldstone bosons is 
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and if^ is a 3x3 matrix with a "1" in the (1, 3) entry, inducing a s — > m tran- 
sition. Octet dominance (A/ = |) has been assumed and thus contributions 
from the 27 component of the As = 1 hamiltonian have been neglected. 
The constant \g^\ = 5.1 is fit to the amplitude for K — * 7ctt{I = 0). 

In computing observables in xPT, the external momentum and quark 
masses are expansion parameters in which the form factors are expanded, 
e.g. G = G'(i) + + Q{3) ^ ____ Ti^g fo^jn factor G'(^) is associated with 
a contribution of order Q^*""^, where Q = p,m, the external momenta or 
meson mass. The same expansion and notation is used for the F± form 
factors. Unlike the contributions from the K2 component, contributions 
from the Ki component are suppressed by a factor of e. However, the 
leading order contribution to — > 7r"'"'7r~7*, r = 1, is from tree-graphs 
involving the Ki component, as shown in fig. (El). A simple calculation 
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Figure 3: T/ie leading final-state interactions in Ki 7r+7r^7*. ^4 so/zc? 
square denotes a weak interaction and a solid circle denotes a strong inter- 
action. These contribution are proportional to e. 



yields 

^^■^ - q' + 2q-p^ ' ^-'^ - q- + 2q.p_ 

and G^^-* = 0, which has the correct symmetry under p± p^ as dis- 
cussed previously. The subscript on the F± form factors indicate that the 
contribution comes from the Ki component. As all constants appearing in 
eq. are determined by other processes, this is a parameter free leading 
order prediction. Final state strong interactions that contribute to F±^i are 
important for CP- violating asymmetries such as -Bcp- The leading final 
state interactions associated with F±^i are generated by graphs shown in 
fig. (0). Retaining only the imaginary parts of the graphs, naively enhanced 
by factors of vr over the real parts, we have 



Im 



p{2) 



^ - ml) (4m|- - 2ml) L 4m| 



+ 2g ■ p+ y mj^ 



which is the leading term in building up e*''°, where So is the I = J = nn 
phase shift evaluated at s = mj^. 

Decay of the K2 component is described by both the G and F±^2 form 
factors starting at r = 2, as can be seen from eq. (|^). At this order in xPT, 
is a constant that must be determined from data. The Ml fraction of 
the decay rate for Ki 7r+7r~7 is reproduced if G*-^-' = 39.3, where higher 
order (momentum dependent) contributions have been neglected, and C-^-* 
is real. The dalitz plot for 7r+7r~7 indicates that there is non- 

negligible momentum dependence in G, and therefore higher order terms 
will be important [in, ITTi. This introduces an uncertainty into the prediction 
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Figure 4: Contributions to 7r+7r^7* from the K and tt charge radii. 

The solid square denotes a weak interaction and the lightly shaded circle 
denotes the sum of one-loop graphs and counterterms that form the charge 
radius of either the K or the tt. 



of differential rates and CP- violating asymmetries at the order to which we 
are working. At r = 3 there are contributions, not only from loop diagrams, 
but also from higher order weak interactions and the Wess-Zumino term. 
However, as before we are able to compute the leading contribution to the 
imaginary part of G, that go to build up the final state interactions, e 



iSi 



where 6i is the phase shift for vrvr scattering in the I 
is found that 



Im 



^(3)1 = ^(2) 
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J = 1 channel. It 



(7) 



where s is the invariant mass of the vr+vr^ system. 

The F-j- 2 form factors do not arise only from the charge radius of the 
as was assumed in the analyses of 0. In fact, the charge radius 
is one of several different types of one-loop graphs arising at r = 2 that 
give rise to g^-dependence in the -F±,2- The diagrams giving contributions 
from the charge radii of the and the vr^ are shown in fig. @. The sum 
of the one-loop diagrams contributing to the charge radius is finite, 
while those contributing to the vr charge radius are divergent and require 
the counterterm 



£ = —i Xc 



(8) 
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Figure 5: One-loop contributions to K2 7r+7r^7* from diagrams that do 
not contribute to the charge radius of the K or tt. A solid square denotes 
a weak interaction and a solid circle denotes a strong interaction. 



where Q is the hght-quark electromagnetic charge matrix, and F^"^ is the 
electromagnetic field strength tensor. The coefficient Acr = —0.91 ± 0.06 
has been determined from measurements of the vr charge radius. 

Diagrams that are not charge radius type contributions are shown in 
fig. (^. Analytic expressions for the diagrams shown in fig. (^), given in 
[§, 1^, are somewhat lengthy and we do not present them here. The sum of 
the graphs in fig. is not finite and the counterterms that enter at this 
order are described by the lagrange density ||12| 



Gf 31 e f^ 
16 V2 



+a2 F^' Tr Q{J:D^T))H^{J:D,T) 
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+a3 F'^'' Tr [iJ^[g, S]D^StSD,Et - i/^D^SD.StSpt, Q] 
+a4 F'^'^ Tr [/7^SD^S"^[g, S]D^S"^j 1 + /i.e. , 



(9) 



where the constants 01,2,3,4 must be determined from data. The combination 

7r+7r^7* is 



of counterterms that contributes to 



w = as - 04 + -(ai + 202) + Acr , 




while the combination that contributes to K 



TT^e^e IS 



2 1 

w+ = -(ai + 202) - 4Acr - - log 

6 D 



1 
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(10) 
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One has the choice to write the F±^2 in terms of w, or to use the known 
values of Acr and oi + 2a2 and define the finite, //-independent combination 
w^L = 03 — ci^^- The value of wl can be determined from the rate for 

The differential decay rate is the incoherent sum of the rates from the 
three form factors, 



dT 



dT 



G 



dq^ 



+ 



dVf 



dq^ 



+ 
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F2 



dq^ 



(12) 



due to the symmetry properties of the amplitudes. In fig. (^) we have shown 

where y = {rriK — 2^,,^), for 



1 dT 
dy ' 



the differential branching fraction 

different values of w^, given the central value of w+ = 0.89||13| and the 
central value of A^ = —0.91. The contribution to the differential rate 
from F-j- 2 vanishes as q^ — >■ 0, but clearly dominates the high q^ region 
(for most values of wl). Except for the ^ region, the contribution 
from G dominates over the contribution from -F±,i. It is clear that in 
order to determine wl a relatively high cut on the e^e~ invariant mass 
must be made. To emphasize this point, the branching fraction for — > 
7r"'"7r~e"'"e~ with a cut of g^^^. > (2 MeV)^ is (using the parameter values 
already discussed) 



Br 



(I6.I + 10.7 + [3.7 - ?,.hwL + 0. 



10" 



(13) 



where the first contribution is from G, the second is from Fi and the third is 
from F2. In contrast, the branching fraction with a cut of q^^^ > (80 MeV)^ 



IS 



Br= (0.6O + 0.07 + [ 1.9 - 1.8 wl + 0.4 wl 



10" 



(14) 



With the presently available branching fraction of Br = (3.32±0.14±0.28) x 
10"^ from KTeV§, which has a g^^^ > (2 MeV)^ cut, wl = 4.7 ± 0.7 or 
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Figure 6: The branching fraction verses y, where y = \^/{mK — 

2m7r). The dot- dashed, dashed and dotted curves are the contributions from 
F±,i, G and F±^2 respectively, while the solid curve is the sum of the contri- 
butions. The three different plots correspond to the counterterm wl taking 
the values 0.40, 0.47 and 0.50 respectively. 



—0.6 ± 0.7, but these values depend sensitively upon G*^^) and F±^i for 
obvious reasons. Only an analysis of the entire differential spectrum, or 
the shape of the tt+tt^ invariant mass distribution will place more stringent 
bounds on Wl- 

One of the most exciting aspects of Tr+vr^e^e^ is the large value 

of Bcp that is predicted |^, ||, || and also recently observed by KTeV0. 
-Bcp is defined to be 



-Bcp = {Sign [sin cos ( 



{Sign 



P++P- 



) 



(15) 



V \p+ + p- 

where is the Pais-Trieman variable depicted in fig. (|^, Ug is the normal 
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to the plane formed by the momenta of the e"'"e~ pair and is the normal 
to the plane formed by the momenta of the tt+tt^ pairs. It is integrated 
over the momenta of the final state particles with any specified cuts. The 
integrand that contributes to Bqp is proportional to the combination 



+ (Fl'l-F^'IY + ...1 . (16) 



The contribution from Re 
does not constitute a comp 



G(3) 



has not been computed, and therefore this 
ete computation of Bqp to next-to- leading order. 
However, the omitted contribution is expected to be small|]^, |^. With a 
cut of q^^^^ > (2 MeV)^ this asymmetry is found to bep|, 3 



Bcp = 9.2% + 4.2 % + ... = 13.4 % , (17) 

with an uncertainty estimated to be of order ~ 2% based on the differ- 
ence between the leading and next-to-leading order contributions. This 
is in complete agreement with the recent KTeV0 observation of Bqp = 
13.6 ± 2.5 ± 1.2 % for this invariant mass cut, and consistent with the cal- 
culations of 1^. The next-to-leading order contribution of 4.3% is from 
the final-state interactions associated with F±,i. It is important to note 
that F-j- 2 does not contribute to Bcp, and hence the uncertainty in deter- 



mining wl does not impact this discussion. As emphasized by SehgallM, 
good agreement between theory and the current experimental value of Bqp 
is obtained within the context of the standard model, with CP-violation 
from e and CPT-conservation. Recent discussions of the implication of this 
observation for T- violating interactions can be found in [^, |l6l. While 
reversing the momenta of the final state particles does change the sign of 
Bqp (it is T-odd), the initial and final states in the decay have not been 
interchanged. Therefore, a direct connection to T-violating interactions is 
absent. 

In conclusion, I have presented a systematic analysis of the decay 
TT'^TT~e~^e~ in chiral perturbation up to next-to-leading order. This analysis 



9 



differs from that of |^ in the form of the 7r+7r~7* dependence upon 

q^. The size of this contribution is determined by a counterterm, wl-, that 
presently is only loosely constrained, but could be determined from the 
existing KTeV data with appropriate kinematic cuts. The large value of 
the CP-asymmetry, -Bcp, that was predicted to arise naturally from e has 
been confirmed by the KTeV collaborationP|. 

I would like to thank Jon Rosner and Bruce Winstein for putting to- 
gether a very stimulating workshop. This work is supported in part by 
Department of Energy Grant DE-FG03-97ER41014. 
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